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Obesity is a primary risk factor for numerous metabolic diseases including metabolic syndrome, type II 
diabetes (T2DM), cardiovascular disease and cancer. Although classically viewed as a storage organ, the field 
of white adipose tissue biology is expanding to include the consideration of the tissue as an endocrine organ 
and major contributor to overall metabolism. Given its role in energy production, the mitochondrion has 
long been a focus of study in metabolic dysfunction and a link between the organelle and white adipose 
tissue function is likely. Herein, we present a novel method for assessing mitochondrial bioenergetics from 
whole white adipose tissue. This method requires minimal manipulation of tissue, and eliminates the need 
for cell isolation and culture. Additionally, this method overcomes some of the limitations to working with 
transformed and / or isolated primary cells and allows for results to be obtained more expediently. In addition 
to the novel method, we present a comprehensive statistical analysis of bioenergetic data as well as guidelines 
for outlier analysis. 

© 2014 The Authors. Published by Elsevier B.V. 
This is an open access article under the CC BY-NC-SA license 
(http://creativecommons.Org/licenses/by-nc-sa/3.0/). 



Introduction 

Metabolic syndrome is a growing epidemic in the United States 
and worldwide. Currently, approximately 34% of the US population 
are living with this diagnosis [1,2] and it is expected that this number 
will continue to increase. By 2030, it is anticipated that the financial 
burden of healthcare for patients requiring treatment for diabetes, 
heart and stroke will be $48-66 billion annually [3]. One of the hall- 
marks of metabolic syndrome is central obesity, an increase in visceral 
white adipose tissue (WAT) around the midsection. Central obesity 
not only contributes to metabolic syndrome but also increases indi- 
vidual risk for type II diabetes mellitus (T2DM) and cardiovascular 
disease (CVD) [4]. 

WAT has classically been viewed as a storage organ, containing 
metabolic substrates in the form of triglycerides as well as insulation 
and cushioning for the internal organs [5]. A more contemporary per- 
spective that has surfaced over the last decade is one of WAT being an 
endocrine organ and a major factor in metabolism and homeostasis 
[6-9]. For example, in addition to triglyceride storage, WAT secretes 
numerous hormones and cytokines (adiponectin, leptin, TNF-oc, etc.) 
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[9]. Coupled with this viewpoint of WAT is the hypothesis that numer- 
ous cardiometabolic pathologies have mitochondrial dysfunction as a 
contributing factor [10-13]. Mitochondria are central to metabolism, 
yet prone to damage that can adversely affect their function [14,15]. 
The assessment of mitochondrial function in WAT has been identified 
as a key target for T2D and obesity [16] and current studies utilize 
both isolated adipocyte mitochondria [17], immortalized adipocyte 
cell lines [18], and isolated primary cells [19]. Isolation of mitochon- 
dria from WAT requires multiple steps and results in measures of 
organelle function in the absence of its native cellular environment. 
Also, functional measures of immortalized/isolated cells that have 
been manipulated to allow for multiple passages may substantially al- 
ter aspects of mitochondrial function. A method for direct assessment 
of mitochondrial function from whole WAT would provide an alterna- 
tive or complementary approach for exploring adipose bioenergetics. 
In this regard, methods for mitochondrial bioenergetic measurement 
from whole adipose tissue have been previously reported [20,21 ] but 
have certain limitations in that certain parameters of mitochondrial 
function could not be defined. 

A methodology herein is provided that utilizes minimal amounts 
of tissue and increased number of timepoint collections to ensure 
tissue equilibration. This approach allowed the measurement of ad- 
ditional parameters of mitochondrial function including ATP-linked, 
maximal oxygen consumption and reserve capacity. Additionally, a 
reproducible data filtering method was employed to account for 
among replicate variance often seen when analyzing whole tissue 

n access article under the CC BY-NC-SA license (http: / / creativecommons.org/licenses/ 



K.J. Dunham-Snary et al. / Redox Biology 2(2014) 656-660 



657 



samples (relative to isolated cells or mitochondria), presumably due 
to differential absorptive properties of whole tissue versus isolated 
cells. Therefore, the utility of the methodology and approach pre- 
sented should aide those interested in evaluating raw bioenergetic 
data obtained from whole tissues including WAT. 

For these studies, mitochondrial bioenergetics was compared be- 
tween two WAT depots, inguinal (iWAT) and epididymal (eWAT) from 
C57BL6/J mice; a mouse model classically used as a control for car- 
diovascular and metabolic disease studies. 

Methods 

Acquisition and preparation of adipose tissue 

12-Week old male C57BL6/J mice were anesthetized by step-wise 
chloral hydrate overdose and euthanized via removal of the heart 
and exsanguination. The Institutional Animal Care and Use Commit- 
tee (IACUC) of the University of Alabama at Birmingham approved 
all procedures involving animals. Inguinal (iWAT) and epididymal 
(eWAT) adipose depots were excised and placed in small petri dishes 
with "rinse media" comprised of DMEM with 25 mM glucose and 
25 mM HEPES. Tissue was gently agitated to remove hair and blood 
and placed in clean rinse media. Each fat pad was sampled six times 
using a 2 mm UniCore Harris Punch (Whatman, Piscataway, NJ), re- 
sulting in 2 mg tissue samples. Attention was made to avoid any visi- 
ble microvasculature during sampling and tissue punches were taken 
from the same section of tissue as close together as possible. Tissue 
punches were kept in fresh rinse media until application to an XF24 
Islet Capture Microplate (Seahorse Bioscience, North Billerica, MA). 

Application of tissue to assay plate 

An XF24 Islet Capture Microplate Screen was loaded onto the 
Screen Applicator Tool and the tool was positioned with the screen 
facing up. A total of 4 mg of tissue (2 tissue punches) from each fat 
pad were placed onto the screen with forceps, resulting in three tissue 
replicates per fat pad. The screen was then snapped into the appropri- 
ate well of the Islet Microplate. Rinse media was added to each well 
after application of the tissue and screen. After applying all tissue 
samples and rinse media, each well was rinsed two additional times. 
After removal of the final rinse, running media (DMEM with 25 mM 
glucose and without HEPES) was added and removed to remove any 
remaining rinse media. Finally, 450 \iL of running media was added 
to all sample and control wells. 

Mitochondrial effectors 

After a series of titration experiments for each mitochondrial ef- 
fector, the following working concentrations were determined to 
elicit maximal responses: 8 \ig/ml oligomycin, 8 \iM FCCP (carbonyl 
cyanide p-trifluoromethoxyphenylhydrazone), 12 \xM antimycin A 
(AA), and 3 \iM rotenone (Rote). The effects of these mitochondrial in- 
hibitors on oxygen consumption rate (OCR - pmol oxygen per minute) 
were assessed approximately 1 h after commencement of basal oxy- 
gen consumption measurements. After multiple test experiments, the 
number of basal measurements was increased to both allow for tissue 
settling and equilibration on the islet capture microplate. 

XF24 Extracellular Flux Analyzer raw data output was transformed 
to OCR using the "Level (Direct) Akos" algorithm in a middle point 
mode, as previously described [22]. Cellular bioenergetic parameters 
including ATP- and non-ATP- linked oxygen consumption as well as 
maximal oxygen consumption and reserve capacity were calculated 
as previously described [23]. Briefly, ATP-linked oxygen consumption 
was calculated as the difference between basal oxygen consumption 
and OCR measured after the addition of oligomycin to inhibit complex 



V (ATP synthase). Non-ATP- linked oxygen consumption was the re- 
mainder of the basal OCR (minus the AA/Rote minimum rate which 
is defined as non-mitochondrial oxygen utilization). Maximal OCR 
was determined by subtracting the AA/Rote minimum reading from 
the maximal reading following addition of FCCP. Reserve capacity 
(or, spare respiratory capacity) is calculated by subtracting the basal 
oxygen consumption rate (respiration rate prior to addition of com- 
pounds) from the maximal oxygen consumption after the addition 
of FCCP (potential maximum respiratory capacity). It is hypothesized 
that reserve capacity exists as a compensation mechanism for cells 
when presented with increased workload or as a stress response [23]. 
A sample OCR trace from WAT is shown in Fig. 1. 

Data summarization and exploratory observations 

Raw OCR trace files were grouped according to eWAT and iWAT, 
and visualized using Microsoft Excel 201 1 . Trace data from each sam- 
ple was divided into five time periods, from which seven OCR vari- 
ables were calculated (basal, mitochondrial, non-mitochondrial, ATP 
linked, non-ATP linked, maximum, and reserve capacity) (Fig. 1 ). Each 
OCR measurement consisted of 3 min of mixing; 2 min wait time and 
3 min of continuous measuring of 0 2 levels. OCR values correspond- 
ing to the equilibration period were ignored, due to physiological 
acclimation of the tissues to the well microenvironment. OCR values 
and slopes for each time period were summarized as the mean value 
across measurements, and the coefficient of least-squares regression, 
in order to assess the distribution of variance among samples. 

Quality control 

Exploratory observations of reserve capacity revealed a high de- 
gree of variance among OCR values within eWAT and iWAT data 
sets (Table 1 ). Reserve capacity was selected as the first variable to 
analyze as a result of the maximal oxygen consumption being the 
most discriminatory during analysis of experiment success. To test 
whether this variance represented biological (among sample) vari- 
ance, or methodological imprecision (among replicates), single fac- 
tor ANOVAs were conducted on both groups using individual mouse 
identification number as the explanatory categorical variable. For this 
and all subsequent statistical analyses, a significance level {a) of 0.05 
was assumed. The ANOVA failed to detect a significant individual 
(mouse) effect on reserve capacity in the eWAT (p > 0.4) and iWAT 
(p > 0.4) groups (Fig. 2). The effect of well position on the XF-24 plate 
was also evaluated as a possible source of measurement bias, and an 
ANOVA was conducted on pooled trace data using well number as 
the explanatory variable. This analysis also failed to detect a signifi- 
cant effect on reserve capacity (p > 0.8). Together, the ANOVA results 
suggested that a majority of among-replicate variance in basal slope 
could not be explained by the well position on the XF-24 plate or 
variation among animals. Therefore it was surmised that the majority 
of observed variance is the product of some unanticipated source of 
methodological imprecision, which is a potential source of noise that 
may impede detection of real biological patterns. 

Consequently, it was necessary to filter data according to a classifi- 
cation method that would differentiate acceptable and unacceptable 
samples. Replicate traces for each mouse were ranked by reserve 
capacity, and pooled by rank (e.g. "maximum", "median", and "min- 
imum"). Reserve capacity values were pooled according to rank to 
compare distribution parameters (symmetry, kurtosis, skew, etc.). 
Exploratory analyses identified a heavy tailed distribution and un- 
acceptable coefficient of variation for pooled reserve capacity data, 
which was likely due to bias in the minimum group (Table 2). This 
suggested that a majority of samples with the lowest reserve capac- 
ity measurements were affected by nonrandom effects, interpreted 
here as interference from unidentified non-biological factors. Right 
skewed data presents a challenge for certain classification schemes, 
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Fig. 1. Example trace from a single XF24 well depicting oxygen consumption rate of C57 iWAT. Bioenergetic parameters including basal, ATP-linked and maximal oxygen 
consumption as well as reserve capacity are outlined on the trace with corresponding mitochondrial effectors. 



Table 1 

Raw OCR data from eWAT and iWAT. 



eWAT 



iWAT 



Pooled 



N 


13 


13 


26 


Basal OCR 


168.66 ± 47.56 


200.27 ± 102.32 


184.46 ± 79.82 


Maximum 


287.28 ± 97.25 


333.32 ± 178.42 


310.3 ± 142.73 


Non-mitochondrial 


77.67 ± 32.45 


103.53 ± 51 


90.6 ± 43.91 


Mitochondrial 


209.6 ± 92.57 


229.8 ± 155.74 


219.7 ± 125.94 


Non-ATP-linked 


125.92 ± 44.94 


158.84 ± 98.32 


142.38 ± 76.76 


ATP-linked 


42.73 ± 35.28 


41.43 ± 28.69 


42.08 ±31.51 


Reserve capacity 


118.62 ± 80.88 


133.05 ± 90.24 


125.84 ± 84.28 



Description: Table summarizing comparison of mean OCR values between raw eWAT and iWAT composite trace data sets. Data represent mean ± standard deviation. All values 
in pmol/min. 
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Fig. 2. Reserve capacity as measured for replicate eWAT (hollow markers) and iWAT 
(filled markers) samples, grouped by mouse. No significant differences were found 
between individual mice. Samples within shaded area were omitted from analyses, 
following the head/tail breaks classification method of Jiang [25]. 



such as Jenk's natural breaks, which seeks to decrease variance within 
replicates and maximize the variance between replicates, but is less 
sensitive to the hierarchy underlying heavy tailed data [24]. 

In contrast, a relatively new technique, head/tail breaks, is suit- 
able for data with this distribution because it emphasizes differences 



among less frequent observations. The algorithm is described in detail 
by Jiang [25], but the concept is as follows. The mean value among 
the data set is used to partition the data into two parts (the "head" 
and "tail"), and the procedure is iterated until the head values are no 
longer heavy tailed. In this case, no further iterations were conducted 
after the mean and median of the head group differed by less than 1 0%. 
The second iteration effectively reduced right skew in the minimum 
and median rank groups, providing a reserve capacity cutoff value 
(58.11 pmol/min). Samples falling below this value were suspected 
to have been influenced by methodological imprecision and omitted 
from further analysis. Composite OCR traces (eWAT and iWAT) for 
each mouse was calculated using values averaged across the filtered 
data set. Composite traces were used in comparative analyses. 

Comparison of OCR between fat pads and with other tissue types 

When pooled across both fat pads, the mean reserve capacity 
(FCCP-basal) (204.96 ± 120.34 pmol/min) was lower than values 
reported for single muscle fibers [26], but similar to values re- 
ported for isolated brown adipocytes [27]. A two-tailed t-test as- 
suming equal variance did not detect a significant difference be- 
tween eWAT( 204.9 ± 102.00 pmol/min) and iWAT reserve capacity 
(205.02 ± 142.89 pmol/min, p > 0.89) (Table 3). 

In addition to these comparisons, it is important to note that 
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Table 2 



Relative variance and normality of raw and filtered data sets. 




N 


Cv 


PPCC 


Raw respiration data 








Maximum 


26 


0.61 


0.97 


Median 


26 


0.84 


0.97 


Minimum 


26 


1.84 


0.96 


Filtered respiration data 








Maximum 


24 


0.55 


0.96 


Median 


16 


0.47 


0.96 


Minimum 


9 


0.47 


0.94 


Description: Table summarizing relative variance and normality compared between raw and filtered trace data setsCv - 


coefficient of variation and PPCC - normal probability plot 


correlation coefficient. 








Table 3 








Filtered OCR data from eWAT and iWAT. 










eWAT 


iWAT 


Pooled 


N 


12 


11 


23 


Basal OCR 


190.97 ± 57.65 


221.43 ± 97.02 


205.54 ± 78.63 


Maximum 


395.86 ± 131.33 


426.45 ± 200.23 


410.49 ± 164.59 


Non-mitochondrial 


75.28 ± 39.91 


120.48 ± 57.31 


96.9 ± 53.13 


Mitochondrial 


320.58 ± 122.40 


305.97 ± 172.23 


313.59 ± 145.02 


Non-ATP-linked 


145.1 ± 59.05 


186.54 ± 85.01 


164.92 ± 74 


ATP-linked 


45.87 ± 42.55 


34.88 ± 24.22 


40.62 ± 34.69 


Reserve capacity 


204.9 ± 102.00 


205.02 ± 142.89 


204.96 ± 120.34 



Description: Table summarizing comparison of mean OCR values between filtered eWAT and iWAT composite trace data sets. Data represent mean ± standard deviation. All values 
in pmol/min. 



looking at WAT in a disease state presents the unique challenge of 
macrophage infiltration wherein a large portion of adipose tissue 
may be macrophages [28]. Under these conditions, it is likely that 
OCR values will be significantly different than those where there is 
not an abundance of macrophages. Mounier et al. report OCRs that 
are orders of magnitude higher than those reported herein from WAT 
[29], even after those results were normalized to protein levels. This 
suggests that during initial screening of data, it will be evident that 
the dominant cell population in any particular well is not WAT, but 
is in fact inflammatory cells. Under these circumstances, it may be 
necessary to normalize bioenergetic data not only to an adipose tis- 
sue marker, but also a macrophage marker. Additionally, histological/ 
immunohistochemical and gene expression analyses may also be war- 
ranted. 

Discussion and conclusions 

With obesity and metabolic disease incidence on the rise world- 
wide, a greater understanding of adipose tissue energetics and its 
contributions to whole body metabolism is needed. The goal of this 
study was to develop a reliable method for assessing mitochondrial 
bioenergetics from whole white adipose tissue by obtaining the great- 
est amount of information with a minimal amount of tissue. 

Herein, we report the development of a method that allows the 
assessment of oxygen utilization in both visceral and subcutaneous 
adipose tissue in mice. In these analyses comparing bioenergetic pro- 
files between fat pad tissues derived from C57BL6/J mice, we found 
no significant differences between tissues obtained from epididymal 
(visceral) and inguinal (subcutaneous) sources of fat. This result how- 
ever is not surprising, as the animals used in this study were fed 
a standard chow diet and were considered lean at the time of tis- 
sue harvest. Since the animals did not have significant levels of WAT 
deposition and were not considered obese, the lack of variation be- 
tween subcutaneous and visceral adipose tissue was not unexpected. 
It is possible that introducing additional variables such as a change in 
diet would alter the bioenergetics of WAT. 

In circumstances where the adipose tissue is from a disease state, 
such as obesity, where the adipocytes are likely to be hypertrophic 



[30,31], an additional normalization would be required. Since larger 
cells (and therefore, fewer mitochondria) would comprise the same 
tissue explant, it is reasonable to assume that bioenergetics would 
be altered. In these cases, one could normalize all bioenergetic data 
to relative mtDNA copy number using a method described by West- 
brook et al. [32]. The data described herein was normalized to relative 
mtDNA copy number, though copy number was not different between 
iWAT and eWAT of lean C57BL6/J mice (data not shown). 

Additionally, the results presented herein represent one strain of 
laboratory mouse and consequently, are limited in scope. However, 
the primary aim of this study was to establish a method for assess- 
ing bioenergetics, and the method development was successful. The 
method presented herein is specific to iWAT and eWAT of lean, chow- 
fed C57BL6/J mice and should be interpreted as specific to both the 
strain of mouse and fat pad. All response rates were determined dur- 
ing titration and time-course experiments for each effector and the 
optimal number of measurements were determined and utilized in 
the final described method. These titration and time-course experi- 
ments (data not shown) allowed for stabilization of OCR after each 
effector was added. Once stabilization was observed, the time of sta- 
bilization was noted and the final method was shortened so experi- 
ments were not unnecessarily long. A lag is present after injection of 
oligomycin (see Fig. 1) to account for variation in stabilization time 
not observed with FCCP or Antimycin A/Rotenone. The authors rec- 
ommend that both titration and time-course experiments for each 
effector be performed when using the method herein on any tissue 
or animal model other than iWAT/eWAT from lean C57BL6/J mice. 

This method is distinct from previously published methodologies 
for whole adipose tissue bioenergetics. These changes included: (1) 
reducing the amount of tissue required per well from 10 to 2 mg in 
order to reduce sample consumption as well as increase the likeli- 
hood of saturating tissue with effectors after injection; (2) increasing 
the number of DMEM rinses prior to loading plate into XF24 to en- 
sure removal of all free lipid released during application of tissue to 
microplate; (3) increasing the number of timepoints for both basal 
as well as post-injection measurements to allow for maximum tissue 
equilibration as well as permeation of effectors into tissue; and (4) 
increasing the concentration of all mitochondrial effectors used both 
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to acquire an oligomycin-induced response previously unreported as 
well as to ensure robust response in response to other effectors. (5) 
Employing a reproducible data filtering method to account for among 
replicate variance in respiration rate. 

We thought a more comprehensive and regimented approach for 
implementing data quality control prior to comparative analyses was 
necessary for these studies due to the anticipated degree of impreci- 
sion observed for whole tissue samples, likely reflective of the differ- 
ential absorptive properties of intact adipose tissue plugs compared to 
those reported for studies using isolated cells or mitochondria. As this 
inconsistency could obfuscate biological differences among individ- 
uals and tissue types, a novel filtering method was used to eliminate 
unacceptable runs from the data set. In this study, reserve capacity 
was selected as the sample exclusion criterion, but the method could 
be applied to any measurement that exhibits unacceptable variance 
among replicates (e.g. basal OCR, ATP-linked OCR, etc.). All data filter- 
ing methods suffer from some degree of subjectivity or bias, and the 
approach used here is not intended to represent the ideal or universal 
technique. The head/tail breaks technique is conceptually transpar- 
ent, easily reproducible, and particularly well suited to right skewed 
data. Therefore, respirometry researchers who obtain raw data with a 
comparable degree of unexplained variance among replicates should 
consider this filtering method before conducting comparative analy- 
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